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Abstract. - We report magnetization measurements and magnetic resonance data for SrTiC>3 
doped by manganese. We show that the recently reported coexistent spin and dipole glass (multi- 
glass) behaviours are strongly affected by the distribution of Mn ions between the Sr and Ti 
sites. Motivated by this finding we calculate the magnetic interactions between Mn impurities 
of different kinds. Both LSDA+f/ and many-body perturbation theory evidence that magnetic 
and magnetoelectric interactions are mediated by Mn s + ions substituting for Ti. We propose two 
microscopic magnetoelectric coupling mechanisms, which can be involved in all magnetoelectric 
systems based on incipient ferroelectrics. In the first one, the electric field modifies the spin sus- 
ceptibility via spin-strain coupling of Mn s + . The second mechanism concerns Mn pairs coupled 
by the position-dependent exchange interaction. 
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Introduction. - SrTi0 3 (STO) and KTa0 3 (KTO) 
doped by manganese have attracted considerable attention 
exhibiting simultaneous spin and dipole glass behaviours 
with large non-linear magnetoelectric coupling [T]-[3]. Such 
"multiglass" systems extend non-trivially the frame of 
conventional multiferroicity and give new perspective for 
studies of the phenomenon and potential application in 
microelectronic devices. 

Both STO and KTO are special representatives of the 
perovskite family of ABO3 materials. They are incipient 
ferroelectrics (IF), i.e. they remain paraelectric down to 
zero temperature, but exhibit very large dielectric permit- 
tivity (~ 20000 and 5000 respectively) at low temperature 
due to the softening of a transverse optical mode that cor- 
responds to B sub-lattice oscillations with respect to the 
almost static rest of the lattice. 

The Mn impurities in STO may substitute both for Sr 
and Ti; they will be denoted as Mn^ and Mils respec- 
tively. Isolated impurities are paramagnetic, Mng be- 
ing an isotropic centre with formal valency Mn 4+ , and 
a spin S = 3/2 [3], while Mn^ which has valency 2+ and 



S = 5/2, is isotropic at T > 100 K, and axial at low tem- 
perature [S|. According to the interpretation of ESR mea- 
surements [3J[5] , which were recently confirmed by density 
functional theory (DFT) calculations 0[7] and EXAFS 
experiments [5], the Mng impurity resides in the octahe- 
drally coordinated cubic position B of the perovskite lat- 
tice, and Mn^4 is displaced from A position thus forming 
electric dipoles in addition to magnetic ones. 

Ceramic samples of STO doped by 2% of manganese 
exhibit spin- and polar- glass properties at temperature 
below T g w 38 K. Moreover, a substantial non-linear mag- 
netoelectric coupling was measured [THS). A similar be- 
haviour for KTO:Mn system was also found [3]. The in- 
teraction of electric dipoles formed by off-central Mn^ im- 
purities has the same nature as the interaction of other 
dipole impurities in IF, its mechanism is rather well un- 
derstood [9]. In this paper we concentrate on magnetic 
and magnetoelectric interactions in STO:Mn. We show 
that the presence of off-central Mn^ + ions substituting for 
Sr is necessary but not sufficient to induce the multiglass 
behaviour, and that the magnetic interactions are medi- 
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ated by Mn B + ions substituting for Ti. 

Some aspects of the considered problem are interesting 
from the fundamental point of view. The interaction be- 
tween ions with different d-shell filling and different spins 
connected by several bridging ligands requires a general- 
ization of superexchange theory. The dependence of spin- 
Hamiltonian parameters on the external electric field is a 
non-trivial application of the ligand field theory for ions 
in a highly polarizable medium. 

Experiment. — For our experimental studies we have 
used two groups of STO:Mn ceramic samples prepared in 
different labs. We will call them I and II. Type-I ceramics 
with the formal chemical composition Sro.9sMno.o2Ti03 
were prepared by mixed oxide technology described else- 
where [10] , In particular, reagent grade SrCC^, TiC>2 and 
MnC>2 were mixed in appropriate amounts, ball milled, 
dried and calcined at 1100 °C for 2 h. The calcined pow- 
ders were again milled, pressed isostatically and sintered at 
1500 °C for 5 h. For the type-II ceramics (formal chemical 
composition Sro.96Mno.o4Ti03), instead of M11O2, powder 
of MnC03 was used. The calcination was performed at 
1150 °C for 4 h and the sintering at 1360 °C for 2 h. Both 
types of samples possess Mn^ + off-central impurities and 
have similar dielectric properties, but as shown below their 
magnetic responses are strikingly different. 

The magnetic measurements were performed using a 
SQUID magnetometre (MPMS-6S Quantum Design) in 
the temperature range 4.5 - 100 K. The zero-field-cooled 
(ZFC) and field-cooled (FC) susceptibilities were mea- 
sured in an applied field of 100 Oe. The main experimental 
result of this paper is illustrated in FigfT] . While magnetic 
properties of type-I ceramics are similar to that reported 
in Refs. |T][2], name ly they exhibit spin-glass behaviour 
at T < 45 K (Fig. la), the type-II ceramics remain para- 
magnetic down to liquid helium temperature (Fig. lb). 
In addition, for the ceramics I the magnetization exhibits 
a hysteretic behaviour with a finite remanence and co- 
ercivity (Fig. lc). On the contrary, the magnetic loops 
measured on ceramics II (Fig. Id) show no hysteresis and 
remanence. In this case, the magnetization is almost en- 
tirely determined by the paramagnetic contribution. 

In order to study the distribution of Mn ions in the 
lattice and their individual magnetic properties -they are 
the only source of magnetism in our system- we have per- 
formed ESR measurements at 9.2 GHz in the standard 3 
cm wavelength range at temperatures from 4.2 up to 300 
K. An Oxford instrument ESR 900 cryosystem was used. 
In both types of ceramics we have found the spectra of 
Mn^ + and Mn B + ions, which were described in detail in 
our previous paper and other early publications (see, e.g. 
Refs. JMS])- The main difference between the two types of 
ceramics revealed by ESR is the ratio of Mn 2+ /Mn 4+ ion 
concentration, which is 70/30 and 90/10 for types I and 
II respectively, as shown in Fig. l.e-h. The absolute con- 
centration of Mn^ + in both ceramics is about lH-1.3 % as 
the ceramics contain also some amount of MnTiO^ inclu- 
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Figure 1: (Coloronline) Comparison of magnetic properties for 
ceramics I (left) and II (right). a),b) ZFC and FC susceptibili- 
ties as a function of temperature. Magnetization as a function 
of applied field, measured c) in type-I ceramics at 80, 40, 25, 
and 5 K, and d) in type-II ceramics at 40, 25, and 5 K. e),f) 
Experimental (dots) EPR spectra, and the corresponding sim- 
ulations (solid lines). g),h) EPR spectra decomposition which 
shows the percentage of various Mn paramagnetic centres. 



sions reflected in ESR spectra by a broad line (see, Fig. 1 
g,h) and observed previously in electron diffraction spectra 
(Ref. PHI). Although in the type-II ceramics about 70% of 
Mn ions contribute to this MnTiOa fraction, it does not 
markedly influence the magnetic properties of the stud- 
ied samples (Figs, lb, Id) due to relatively small magnetic 
anomaly of MnTiOa at the antiferromagnetic phase tran- 
sition at 63 K (Ref. [TT]). 

DFT calculations. — We have theoretically stud- 
ied the interactions between Mn ions using LSDA+t/ and 
many-body perturbation theory. We have considered the 
Mn impurities when they occupy nearest neighbour posi- 
tions, and looked at the various pairs Mn^-Mn^, Mn^- 
Mng, and Mns-Mns. The Mn^ ions were shifted in var- 
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Table 1: Hamiltonian (1) parameters calculated with LSDA+J7 
for various configurations of Mn^ + -Mn^ + pair in a 2 x 2 x 3 60 
atom supercell (symmetry group P4mm, #99); U =4 eV 

di,A d,2,A E n , meV Jaa, nieV 












-0.008 


0.82 


0.82 


-695.5 


0.23 


0.35 


-0.35 


-26.6 


0.46 


-0.47 


0.47 


-27.2 


0.04 



Table 2: Hamiltonian parameters for Mn 

lated for a 2 x 2 x 2 40 atom supercell (PI, #1); U 



A + - 



Mn^ + pair 









0.5 


-0.64 


-144.6 


1.65 


0.64 


-123.02 


0.51 



calcu- 
4eV 



d A , A E n , meV J AB , meV JabA'b, K 

5.8 
19.1 
5.9 



ious directions from the symmetric A position. 

Density functional theory calculations were performed 
using the full- potential local-orbital (FPLO) code |12j . 
The total energies for different configurations were ob- 
tained, and the results were mapped onto an Mn-Mn pair 
effective Hamiltonian of the form 



F = S„(d 1 ,d 2 ) + J 12 (d 1 ,d 2 )S 1 S 2 , 



(1) 



where d^ is the shift from cubic position. E n is a non- 
magnetic spin-independent interaction, treating the ion 
motion classically, it is a c-number in our approach. The 
spin part of the interaction is of Heisenberg-type. 

We present the results for Mn^"-Mn^ + pairs in Table [T] 
Various Mn^ + displacements along the symmetry axis Z 
were considered: d\ = d% (the ions shift conserving the 
distance between them), d\ = — d% > (ions shift towards 
each other), d\ = — d<i < (away from each other) . In each 
case the results correspond to the total energy minima 
with respect to displacement. The non shifted distance 
between two A sites is 3.9 A [BJ. 

The third column shows that the configuration with par- 
allel shift of Mn ions (di = d 2 = 0.82A) is separated from 
the others by such a large energy, that these will not be ob- 
served within the physically relevant temperature range. 
A similar situation was described in Ref. [TS] for Li-Li pairs 
in KTaO v -j. The last column of Table Q] represents one of 
the main theoretical findings of this work, namely, it shows 
that the magnetic interaction of nearest neighbour Mn^ + 
ions, for the most probable parallel configuration, satis- 
fies |Jaa/&b| < 3 K. In fact this value is on the verge 
of precision of DFT calculations. Despite the large spin 
value of Mn^ + ions, this interaction cannot be responsible 
for magnetic susceptibility anomalies at T ~ 40 K that we 
observe for type I ceramics (FigJTJ) , and that were reported 
in previous studies of STO:Mn [lj[2]. 

Table [5] shows the results for a Mn^ + -Mn^ + pair. Here 
the Mn^ + ion shift was taken to be the same as the one 
found for the isolated ion jBJ; for dA > (<)0 the ions get 
farther (closer) . Finally, for the pair Mn.g -Mng~ we have 
found Jbb = 9.3 meV, i.e. JbbI^-b = 107.9 K. So, we may 
conclude that the presence of Mn^ + ions may be responsi- 
ble for the observed anomalies in magnetic susceptibilities. 

Theory of superexchange. — To reach a better 
understanding of the exchange mechanism we have per- 



formed analytic calculations of the superexchange interac- 
tion between Mn ions within fourth-order many-body per- 
turbation theory. Using resolvent method [TBJ we obtained 
a general formula for the magnetic coupling between two 
spins S± and S2 in terms of hopping integrals U m p n be- 
tween cations (im) (i specifies the cation, m the orbital 
index) and neighbouring ligands (j3n) ((3 specifies the lig- 
and, n the orbital). Taking the hole point of view (no 
fermion on the ligand p orbitals in the ground state) we 
have established Eq.©, where the sum over the cation 
orbitals m and m! is restricted to ground-state occupied 
ones (abbreviation "occ"). A^„ is the energy of the ex- 
cited state (measured with respect to the ground state), 
where one fermion has moved from (im) to (f3n), while Ay- 
is the difference of energy between an excited state with 
TV, — 1 and Nj + 1 fermions on cation i and j respectively, 
and the ground state (GS) with N and Nj fermions per 
cation respectively. If the cations are of the same type, 
one has Ay- = Aji and Aip n does not depend on i. The 
three added terms inside the bracket in Eq.([2]) correspond 
to different paths: the first two involve intermediate ex- 
cited states with (JVj — 1, Nj + 1) fermions on the cations, 
while the third term corresponds to an excited state with 
two fermions on the ligands with repulsion U p when the 
holes meet on the same ligand. 

The parameters U m p n , Aip n and Ay- were extracted 
from the analysis of photoemission experiments reported 
in Refs. [171115] . We denote for the BB-pair: A 12 = Uf fJ , 
A i0n = Af //; for the AA-pair A 12 = Uf ff , A lfSn = Af ff ; 
while for the BA-pair A 12 = A B ff - Af ff + Uf ff , A 2 i = 
A e ff - Af /y + U? f/ , A Wn = Af ff , and A 2fin = Af ff , 
where A*„ = A, + J m [7/9{N t - 1)-Pi] and U l ef} = 
Ui + J m [N + 11/9 - 2j>i], with J m = 5/2Bi + d, N 
is the number of holes in the ground state of Mn.; ion 
and pi is the number of doubly-occupied orbitals. We use 
the following experimental values (in eV) for Mng: Bi = 
0.132, C, = 0.610, Ai = 2 ± 0.5 and U % = 7.5 ± 0.5 (Table 
I of Ref. Q7J,), while for Mn A , B % = 0.119, d = 0.412, 
A, = 7 and U t = 5.5 (Table I, II of Ref. [B]). This 



gives A^ ff 



9.2eV, Af ff 



4.5 ± 0.5eV, Uf ff = 9.9eV 



and U^ff =- 11.6 ± 0.5eV. The hopping parameters are 

expressed in terms of Slater- Koster parameters [TO] V dll j 
for Miiyj^i, they essentially depend on the cation-ligand 
distance, we take V^/V^ = -2.16 [TOJ. Following Ref. 
[T5] , we have assumed that the hopping which links states 
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Table 3: Superexchange values jba = 
Jaa/Jbb with Jbb = (9.5 ± 1.5)(V p B dw ) 



Jba/Jbb, and j'aa = 
4 for different A-shift. 



<i\ 



JBA 



Jba, 
LSDA+U 



JAA 



JAA, 

LSDA+C/ 






0.39 


0.054 


0.008 


-0.001 


0.82 






0.008 


0.025 


-0.64 


0.52 


0.177 






0.64 


0.29 


0.055 







with N + 1 and N holes on some cation, is reduced by 
a factor R compared to hopping linking states with N — 
1 and N holes. The results (and their comparison with 
LSDA+[7) are presented in Table [3] for the parameters 
R = 1.2 and U p = 4 eV, and the shifts are the same as 
those found in LSDA calculations. 

As can be seen, the agreement is qualitative between 
perturbative calculations and DFT results. In addition, 
the super-exchange theory enables to discuss various ten- 
dencies and contributions. The different exchange cou- 
pling Jbb, Jab and J A a are proportional to (V^„.) 4 , 
O^L^O^L) 2 and (^L) 4 respectively A crude estimate 
based on the different distances, leads to Jab/ Jbb ~ 0.1 
and Jaa/Jbb ~ 0.01. One can also expect a predomi- 
nant value for Jbb due to geometry: it corresponds to a 
180° Mn B -0-Mn B link, while the other two are 90° ones. 
The difference between exchange couplings are also due 
to different number of exchange paths and to difference in 
charge transfer values Aj 4 ^ > Af ff . The number of paths 
for J aa is larger than for the other two, since four ligands 
and many orbitals come in, against three ligands for Jab 
and only one for Jbb- Two or more ligands raise the op- 
portunity to have ferromagnetic contributions, which ex- 
plains reduction of Jaa compared to the other two, further 
than what was expected just by distance effect. 



Mechanisms of magnetoelectric coupling. 

We propose two possible mechanisms involving Mn B 
ions, which may be relevant for IF doped by man- 
ganese. In the following we consider the solid solution 
Sri-^Mn^Tii-yMiiyOs. The magnetoelectricity implies 
the dependence of the magnetic susceptibility on an elec- 
tric field. Up to second order, the magnetic susceptibility 



may be written 

d 2 F 
Xlj ~ ~ dHidHj 



H=0 



XO.ij + Xl,ijkEk + X2,ijklEkE[- 

(3) 

where the free energy density F, the magnetic and electric 
field components, Hi and Ei are measured in CGS units. 
To translate this equation in terms of Ref. [5] notations 
with SI units, one should use the relations 



Pijk = (4tt) 2 no\/eoX\,ijk 

Stjki = 8ir 2 no£oX2,ijki ~ 8.78 • 10~ lo X2,ijki sm/VA(5) 



1.67-10- 10 x lj y fe s/A(4) 

16, 



The first mechanism is a one-spin effect: The polariza- 
tion P of the lattice, is accompanied by a lattice strain, 
which is proportional to the square of polarization |20| . 
The Mn B + ion has a 3d 3 configuration and a 4 F GS, which 
is split by the ligand field in contrast to the 3d 5 configu- 
ration and 6 S GS for Mn^ + . In a cubic field the Mn B + ion 
has a fourfold degenerate 4 A 2 GS. When the local symme- 
try becomes axial, an additional splitting arises. Its mag- 
nitude depends on the polarization via the strain. This af- 
fects the magnetic susceptibility. In a paraelectric phase, 
the changes are proportional to the square of the external 
electric field, but in the presence of a net polarization P r 
in a polar phase, a linear dependence appears. Let us note 
that this mechanism will be effective in any ferroelectric 
perovskite doped by paramagnetic ions located at B sites 
and having the 3d 3 configuration. 

The contribution of this mechanism to the magnetoelec- 
tric susceptibility is given below by Eqs.fJSJ) and (fTU)) . Sub- 
stitution of numerical values relevant to the experimental 
conditions of Ref. [2] gives 

X i, zzz « -2.2-10- 9 , /3«-3.6-10- 19 s/A, (6) 
X2, zzzz « -1.2 -10- 10 , c5«-0.1-Kr 24 sm/VA.(7) 

In these estimates, we have used x + y = 0.02, and 
x/y « 70/30, as found in our ESR experiment. The re- 
ported experimental values are [2] f3 f=a —3 • 10 -19 s/A and 
S « —9 • 10 -24 sm/VA. We see that the paramagnetoelec- 
tric susceptibility j3 is in fairly good agreement, but the 
biquadratic coefficient is underestimated by this mecha- 
nism which cannot be the only argument put forward. 

The second mechanism concerns Mn^-Mne pairs. As 
can be seen from Tables [5] and [31 the Mn^ + has non- 
equivalent equilibrium positions in the cell where Mn B 
is present. The superexchange interaction between Miig + 
and Mn^ + strongly depends on the displacement 4a- 
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Compared to the closer situation, the case with ions far 
from each other corresponds to a higher energy AE. As 
shown in detail in the following section, the lattice polar- 
ization increases the number A^-Pr, E) of Mn^ + ions lying 
farther from Miig + ions. This leads to a positive Xi,zzz- In 
the next section we show that this mechanism seems to be 
unimportant for the system studied in Ref. [2], but it may 
be very effective for a system where the interaction energy 
between electric Mn^ + -dipole and P r is comparable with 
AE. 

Calculation of non-linear magneto-electric sus- 
ceptibility. — Here we quantitatively consider the ideas 
outlined in the previous section. This section is a little 
technical and can be safely skipped by someone not inter- 
ested in the details of computation. 

The first mechanism. In order to estimate its contri- 
bution, we follow the ideas of Ref. [21]. The Mn B + ion is 
much more sensitive to the local strain than Mn^ + . In an 
octahedral coordination it has a 4 A2 orbital singlet ground 
state, and the effective spin 3/2 Hamiltonian has the form 



Hr 



where g^ v 



^B9z Z HzS Z + ^Bg±(H X S X + HySy) 



D 



s 2 z -\s{s + i) 



— Ms 2_^ ^iivH,j,H v ,(8) 



/" y 



v <^o|L^|»)(n|L 



2A/cA /u/ , D 

'I'M 



-A 2 



A, 



A, 



\+*-ZZ J1 M/j 1 *-[IV 

g s = 2.0023 is the spin gyromag- 
netic ratio, A the spin-orbit coupling constant, IV'oHI 7 '*)) is 
the ground(excited) state, Eq and E n are the correspond- 
ing energies, I/ 1 is the /i-th component of the orbital mo- 
ment operator, k is the orbital reduction factor. Using the 
ligand field theory |22ll23] , which perturbatively takes into 
account p — d hybridization between paramagnetic ion and 
surrounding ligands, we may obtain |24| 

A zz = -Ag cub (1 - Ke 33 ) /2A, 



A, 



-Ag cub (1 - 5«;e 33 /2) /2A, 



where Ag cub = g cub — g s refers to the g -value for undis- 
torted cubic lattice, k ~ —3.5, —4 is the exponent of p — d 
hopping dependence on distance |19j . The dependence on 
electric field of the previous equations is only in the strain 
633. For multiglass samples we may assume that both, the 
net polarization P r and the electric field are directed along 
the axis Z, while P = P r + (e — 1) E/47T. With an exter- 
nal magnetic held also applied along the Z-axis, we com- 
pute the partition function Zb and the free energy density 
Fb = —(y/v c )0\nZB, with y/v c the density of Mn B + ions, 
v c the host lattice unit cell volume and 9 = ksT. From 
this free energy, we finally derive the linear contribution 
to the magnetic susceptibility (Eq.Q): 



Xi, 



.'/ 



H%K, 



A 



5k 



3g 2 zA 
20 



-2 



(9) 



where we have taken into account that D/9 <C 1 for 
T = 10 K. The (double-)prime indicates the (second-) 



derivative with respect to E. Similarly, for the second 
order contribution, we have 



X2„ 



w% 



^{k(\\' zz f[l0k+l2g zz \/9]/e 

Vr L 



A" 



9zzX Uk + 3gzzX 



(10) 



When we substitute k = —3.5, g cu b = 1.992 |,A« 135 
cm -1 , [25], k « 1, SrTi03 lattice parameter a = 3.9 A, 
(v c = a 3 ), and c/a w 1.002 ©, e(T = 10 K) w 1500, and 
a net polarization P r w 0.7, /zC/m 2 = 2100 esu/cm 2 |10) . 
we obtain the values presented in Eqs. (O and (JT]). 

The second mechanism. It involves a pair of Mn^- 
Miis ions. There are six available positions for Mn^, three 



AE 



Figure 2: Changes of MnA-Mns pair energies in polarized 
medium. Here Ei oc = E + -zrjP is positive, however it can 
be negative, depending on the relative direction of polarization 
with respect to the pair MnA-Mng. The Lorentz factor 7 ac- 
counts for the deviation of the local field from the simple cubic 
case, p m Za&a is the dipole moment of Mn^ + , Za ~ 2 is its 
dynamic charge. 

of them are closer to Mn B and lower in energy: the num- 
ber of Mn^ + in deep (shallow) wells is Nim- Writing 
Ni = m + 7i2 + 113, and N2 = n± + TI5 + n§, one has in 
absence of polarization N2/N1 = exp(—AE/9). However 
with polarization, there is sonre lift of degeneracy, as can 
be seen on Fig. [2J this leads to a redistribution of level 
occupancies, and Nu2) acquires a dependence on the net 
polarization P r and on the external electric field E. With 
Ni + N2 = xyz, ( z = 8 is the coordination number), the 
susceptibility reads 



X = Xo +Xixyz + (X2 -Xi) N 2(P r ,E), 



(11) 



where Xo comes from contribution of everything but Mn^- 
Miib pairs, while x.j{J = 1)2) is the susceptibility of a 
pair Sa = 5/2 — Sb = 3/2 coupled by the exchange in- 
teraction Jab{<1a)i which depends on the Mn^ position 
[2"rJ] . The dependence of x on the electric field comes from 
N 2 {P r ,E). Defining n(t) = (2 + e~*)/(2 + e'), one has 



N 2 (P r ,E) 



xyz 



1 



1 



(12) 

with t = pEioc/9, pEi oc ra cE + A r , c = ps^f/3 and A r = 
4:TrjpP r /3, (e>l was used). See the caption of Fig.[2]for 
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definitions. Finally, making a limited development up to 
second order in electric field E, we obtain 

Xl,zzz = (X2-Xl)Wr,0), (13) 

X2,«« = (Xs-Xi)NS(P r ,0). (14) 

Substituting the values from Table [2J A£ s=s -123 + 
144.6 = 21.6 meV, Ji « 1.7 meV, J 2 « 0.5 meV, and 
7 w —0.2 [9], this gives A r s=s —14 meV, and a very 
tiny contribution of this second mechanism to the sus- 
ceptibilities xi ~ 7.1 • 10~ 15 , %2 ~ 1-6 • 10~ 15 . However, 
since N^ is strongly dependent on AE — A r , the second 
mechanism could be of the same order than the first one 
or even exceed it for AE ps A r (then \i ~ 4.6 • 10 -7 , 
\2 ~ 1.7 ■ 10~ 7 ): this equality can be realized for other 
concentrations of manganese in SrTiC>3 host, or for other 
systems (e.g. Ki_ a; Mri a; Tai_.yMn2,03 [3]). 

Summary. — We have compared the magnetic prop- 
erties of two types of ceramic samples of manganese doped 
SrTiC>3. Based on the data of ESR measurements we 
conclude that the spin-glass behaviour is observed only 
in samples containing an appreciable percentage of Mn B 
ions substituting for Ti, in addition to Mn^ + substitut- 
ing for Sr. Using LSDA+C7 supercell calculation we have 
shown that the exchange interaction between Mn^ im- 
purities is an order of magnitude smaller than those for 
Mnjj + -Mn^ + and Mn^-Mn^" pairs. The analytic many- 
body calculations have shown that the reason for this dif- 
ference is the interference of various exchange paths for 
the Mn^ + -Mn^ + pairs combined with different geometry, 
Mn-0 distances, and stability of the Mn^ + ground state 
configuration in comparison to other pairs. We conclude 
that the presence of Mn B + ions is essential for the forma- 
tion of a collective magnetic state at low temperature. We 
propose two microscopic mechanisms of magnetoelectric- 
ity in SrTiOsiMn which involve Mn fl + ions. 
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with E s = \Jab[S{S + 1) - Sa(Sa + 1) - Sb(S b + 1)], 
the magnetic energy of the pair corresponding to the total 
spin S, \s = in-Bgs) S (S + 1) /(3v c 8) is the correspond- 
ing magnetic susceptibility, and the gyromagnetic ratio is 
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